home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dgtsv.z / dgtsv
Encoding:
Text File  |  2002-10-03  |  3.8 KB  |  133 lines

  1.  
  2.  
  3.  
  4. DDDDGGGGTTTTSSSSVVVV((((3333SSSS))))                                                            DDDDGGGGTTTTSSSSVVVV((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DGTSV - solve the equation  A*X = B,
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE DGTSV( N, NRHS, DL, D, DU, B, LDB, INFO )
  13.  
  14.          INTEGER       INFO, LDB, N, NRHS
  15.  
  16.          DOUBLE        PRECISION B( LDB, * ), D( * ), DL( * ), DU( * )
  17.  
  18. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  19.      These routines are part of the SCSL Scientific Library and can be loaded
  20.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  21.      directs the linker to use the multi-processor version of the library.
  22.  
  23.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  24.      4 bytes (32 bits). Another version of SCSL is available in which integers
  25.      are 8 bytes (64 bits).  This version allows the user access to larger
  26.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  27.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  28.      only one of the two versions; 4-byte integer and 8-byte integer library
  29.      calls cannot be mixed.
  30.  
  31. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  32.      DGTSV solves the equation A*X = B, where A is an n by n tridiagonal
  33.      matrix, by Gaussian elimination with partial pivoting.
  34.  
  35.      Note that the equation  A'*X = B  may be solved by interchanging the
  36.      order of the arguments DU and DL.
  37.  
  38.  
  39. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  40.      N       (input) INTEGER
  41.              The order of the matrix A.  N >= 0.
  42.  
  43.      NRHS    (input) INTEGER
  44.              The number of right hand sides, i.e., the number of columns of
  45.              the matrix B.  NRHS >= 0.
  46.  
  47.      DL      (input/output) DOUBLE PRECISION array, dimension (N-1)
  48.              On entry, DL must contain the (n-1) sub-diagonal elements of A.
  49.  
  50.              On exit, DL is overwritten by the (n-2) elements of the second
  51.              super-diagonal of the upper triangular matrix U from the LU
  52.              factorization of A, in DL(1), ..., DL(n-2).
  53.  
  54.      D       (input/output) DOUBLE PRECISION array, dimension (N)
  55.              On entry, D must contain the diagonal elements of A.
  56.  
  57.              On exit, D is overwritten by the n diagonal elements of U.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDGGGGTTTTSSSSVVVV((((3333SSSS))))                                                            DDDDGGGGTTTTSSSSVVVV((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      DU      (input/output) DOUBLE PRECISION array, dimension (N-1)
  75.              On entry, DU must contain the (n-1) super-diagonal elements of A.
  76.  
  77.              On exit, DU is overwritten by the (n-1) elements of the first
  78.              super-diagonal of U.
  79.  
  80.      B       (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS)
  81.              On entry, the N by NRHS matrix of right hand side matrix B.  On
  82.              exit, if INFO = 0, the N by NRHS solution matrix X.
  83.  
  84.      LDB     (input) INTEGER
  85.              The leading dimension of the array B.  LDB >= max(1,N).
  86.  
  87.      INFO    (output) INTEGER
  88.              = 0: successful exit
  89.              < 0: if INFO = -i, the i-th argument had an illegal value
  90.              > 0: if INFO = i, U(i,i) is exactly zero, and the solution has
  91.              not been computed.  The factorization has not been completed
  92.              unless i = N.
  93.  
  94. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  95.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  96.  
  97.      This man page is available only online.
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.